Kerr-Schild Method and Geodesic Structure in 
Codimension-2 Brane Black Holes 

o 

Bertha Cuadros-Melgar* 

O 

CN Departamento de Ciencias Ffsicas, Universidad Andres Bello, 

Avenida Repiiblica 252, Santiago, Chile 

^ Susana Aguilar^ and Nelson Zamorano* 

, , Departamento de Ffsica, Facultad de Ciencias Ffsicas y Matematicas, 

qh Universidad de Chile, Avenida Blanco Encalada 2003, Santiago, Chile 

u 

^ Abstract 

We consider the black hole solutions of five-dimensional gravity 
CSj with a Gauss-Bonnet term in the bulk and an induced gravity term 

on a 2-brane of codimension-2. Applying the Kerr-Schild method we 
CN derive additional solutions which include charge and angular momen- 

^ turn. To improve our understanding of this geometry we study the 

Q geodesic structure of such spacetimes. 



X 



PACS numbers: ll.25.-w, 04.20.Fy, 04.50.Gh, 04.70.-s 



* e-mail: berthaki@gmail.com 
^e-mail: saguilar@ing.uchile.cl 
■^e-mail: nzamora@dfi.uchilc.cl 



1 



1 Introduction 



One of the unsolved questions of braneworld cosmology is the existence of 
localized black holes on the brane. This puzzle has been investigated al- 
most since the appearance of these alternative-to-classical-gravity models. 
In codimension-1 scenarios the natural first proposal was to consider the 
Schwarzschild metric and study its black string extension into the bulk pp. 
Unfortunately, as intuited by the authors, this string is unstable to classical 
linear perturbations known as Gregory-Laflamme (GL) instability j2]. Since 
then there has been an intensive research to find a full metric by using numer- 
ical techniques [3] or by solving the trace of the projected Einstein equations 
on the brane pi] with special Ansatz [5] or making certain assumptions on 
the projected Weyl term coming from the bulk |6J. A lower dimensional 
version with a Banados-Teitelboim-Zanelli (BTZ) [7] black string was also 
considered in [8]. This solution was obtained from the so-called C-metric. 
The thermodynamical analysis showed that the string remains stable when 
its transverse size is comparable to the 4-dimensional AdS radius and can 
be unbalanced by a GL instability above that scale, breaking up to a BTZ 
brane black hole. 

In codimension-2 the first attempt was proposed in [U] as a generalization 
of the 4-dimensional Aryal, Ford, Vilenkin black hole [TU] pierced by a cos- 
mic string. The rotating version was also presented in and a complete 
study of the gray-body factors was considered in [12]. Another proposal came 
by considering a 5-dimensional bulk with a Gauss-Bonnet term and induced 
gravity on the brane [13] . The solutions are basically of three types. The 
first one is the familiar BTZ black hole which can be extended into the bulk 
with a regular horizon. The second one adds a short distance correction 
and describes a BTZ black hole conformally coupled to a scalar field. There 
is a third solution family that can accommodate any brane metric coeffi- 
cient n(r) provided that it can yield a physically acceptable brane energy- 
momentum tensor. The corresponding generalization to 6-dimensional bulks 
was worked out in [14], where it is shown that the only possible solution, 
namely a Schwarzschild- AdS brane black hole, needs matter in the bulk. For 
all these solutions the Gauss-Bonnet term plays a fundamental role lead- 
ing to a consistency relation that dictates the kind of bulk or brane matter 
necessary to sustain a black hole on the brane. 

The arbitrariness of n(r) in one of the above mentioned solutions is an 
interesting feature that deserves more discussion. This fact motivated us to 
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look for more solutions that could exhibit physically acceptable brane energy- 
momentum tensors. Specifically, we wondered if charged and rotating metrics 
can be included in this general n(r), in particular, taking into account that 
the rotating BTZ black hole has non-diagonal elements that were not included 
in the original derivation of the solutions found in [T3]. With this purpose 
we recur to the so-called Kerr-Schild method. 

The Kerr-Schild coordinates appeared for the first time when obtaining 
the Kerr metric starting from a flat empty space. Later on, A. H. Taub [15] 
generalized this method in order to obtain new solutions by adding to a 
background metric a term proportional to a null geodesic vector and a scalar 
function H. The resulting metric is not a coordinate transformation but 
a new spacetime with a different geometry. The selection of H, although 
arbitrary, must fulfill certain criteria. Not only it has to satisfy the Einstein 
equations when inserted in the new metric but it must also yield a physically 
meaningful energy- momentum tensor. The advantage of this Generalized 
Kerr-Schild (GKS) method is that the new Einstein equations are linear 
in H when written in their covariant and contravariant components. We 
should also mention that the resulting metric can have non-diagonal terms 
that can be cast off by means of a coordinate transformation when no angular 
momentum is involved. The GKS transformation appears as a useful tool in 
astrophysics where, for example, it can generate a singularity free metric that 
describes orbits close to a Kerr black hole event horizon [IB]. Moreover, this 
kind of metric can be applied in numerical relativity to outline the geometry 
at the horizon and even extend it to the interior of the black hole, given 
the absence of singularities of these coordinates in this region. The GKS 
transformation is also used to identify apparent horizons in boosted black 
holes due to its invariance under Lorentz boost [17J. Other solutions were 
also worked out in [18]. This method has also been used to find exact vacuum 
solutions in the context of multidimensional gravity [19|. 

In this paper we applied the GKS method to extra dimensions. The back- 
ground solutions are the black hole metrics described previously. Due to the 
Gauss-Bonnet contribution, the Einstein equations are not linear anymore 
even when written in their covariant-contravariant form, in fact, they are at 
most quadratic in the function H . However, as we will see in this work, the 
method proved to be successful since the equations are still solvable, and it 
was possible to obtain more solutions including charge and angular momen- 
tum. It was also viable to pass from BTZ metrics (pure or charged) to the 
so-called corrected BTZ one, which includes a brane scalar field. 
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Given all these solutions it is pertinent to explore the main properties 
of such spacetimes. In General Relativity one possible way to understand 
the geometrical aspect of the gravitational field is to study geodesies in the 
spacetime permeated by this field. Over the years, the motion of massive 
and massless test particles in background geometries of various higher di- 
mensional theories of gravity has been investigated [20J. In this work we 
analyse the timelike and null geodesic motion in the brane black hole space- 
times mentioned above. With this aim we solve the Euler-Lagrange equation 
for the variational problem associated with the corresponding metrics. The 
set of orbits turned out to be very rich. 

The paper is organized as follows. In section 2 we make a brief review of 
the brane black holes obtained in [13]. Section 3 is devoted to the application 
of the GKS method to these models displaying the corresponding results and 
new solutions. In section 4 we present a complete study of timelike and null 
geodesic behaviour in the whole family of backgrounds. Finally, in section 5 
we discuss our results and conclude. 

2 BTZ String-like Solutions on Codimension- 
2 Braneworlds 

We consider the following gravitational action in five dimensions with a 
Gauss-Bonnet (GB) term in the bulk and an induced three-dimensional cur- 
vature term on the brane [T3] 

+ r\ J Af^) ^ + J <fxC bulk + J <PxC brane , (1) 

where a (> 0) is the GB coupling constant, r 2 = M 3 /M| is the induced 
gravity "cross-over" scale, which marks the transition from 3D to 5D gravity, 
and M 5 , M 3 are the five and three-dimensional Planck masses, respectively. 

The above induced term has been written in the particular coordinate 
system in which the metric is 

ds\ = p)dx^dx v + dp 2 + b 2 (x, p)d6 2 . (2) 

Here g^ u (x, 0) is the brane metric, whereas x^ denotes three dimensions, 
/i = t, r, 0, and p, 9 denote the radial and angular coordinates of the two 
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extra dimensions (p may or may not be compact, and < 9 < 2tt). Capital 
M, N indices will take values in the five-dimensional space. 

The Einstein equations resulting from the variation of the action (HI) are 



(5)N , J2ri(Z)v a n N 

2nb atlM ~ Ml 



Lr M + r c Cr^ g M g v -^—^ - ati M 
where 



rp(B)N . rp(br)v U „N 
± M I ± i, MM 



2nb 



(3) 



ttN 



_L4fi>( 5 ) H NP i4P< 5 ) N j3KP 9 d(5) t,NKLP~\ (a\ 
+ 4 %P rt (5) + 4/ %MP n (5) ~ * n MKLP n (5) \ ■ \V 

To obtain the braneworld equations we expand the metric around the 
brane as 

b(x,p) = P(x)p + 0(p 2 ) . (5) 

At the boundary of the internal two-dimensional space where the 2-brane is 
situated the function b behaves as b'(x, 0) = /3(x), where a prime denotes 
derivative with respect to p. In addition, we demand that the space in the 
vicinity of the conical singularity is regular, i.e., <9 M /3 = and d p g^ u (x,0) = 
[21] • The extrinsic curvature in the particular gauge g pp = 1 that we are 
considering is given by K^ v = g' . Using the fact that the second derivatives 
of the metric contain ^-function singularities at the position of the brane, the 
nature of the singularity gives the following relations [21] 



= — (1 — b') ^— + non — singular terms , (6) 



b 

b 



h non — singular terms . (7) 



From the above singularity expressions and using the Gauss-Codacci equa- 
tions, we can match the singular parts of the Einstein equations (|3| and get 
the "boundary" Einstein equations, 

G$ = j^TW + 2n(l-P)^ gilv . (8) 

We look for black string solutions of the Einstein equations ^ using the 
five-dimensional metric ^ in the form 

dsl = f(p) \-n(r) 2 dt 2 + n(r)- 2 dr 2 + r 2 d<f) 2 } + dp 2 + b 2 (p)d6 2 , (9) 
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where we have supposed the existence of a localized (2+1) black hole on the 
brane, whose metric is given by 



n(r) 2 dt 2 + n(r) 2 dr l + r 



2j„2 , Jljil 



(10) 



In the bulk we consider only a cosmological constant A5. Then, from the 
bulk Einstein equations 



A, 



Ml 



9mn 



By combining the (rr, 00) equations we get 



h + nn 



nn 



1 -4a- 



r J \ b J 

while a combination of the (pp, 99) equations gives 



3-4 



a 

T 2 



( 9 nn „. 9 \ 

(n 2 + nn + 2— + 3f 2 ) 



(11) 



(12) 



(13) 



where a dot implies derivatives with respect to r. The solutions of the equa- 
tions (12) and (13) are summarized in Table [l] [T3] . 
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Table 1: BTZ String-Like Solutions in Five-Dimensional Braneworlds of 
Codimension-2. 

In this table I is the length of three-dimensional AdS space. The BTZ 
solution is given by [7] 

2 

n i(r) = -M + T - . (14) 
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When its mass is positive, the black hole has a horizon at r = l\/M, and the 
radius of curvature of the AdS^ space I = (— A-^)^ 1 ^ 2 provides the necessary 
length scale to define this horizon. For the mass — 1 < M < 0, which is 
dimensionless, the BTZ black hole has a naked conical singularity, while for 
M — — 1 the vacuum AdS 3 space is recovered. 

The corrected BTZ solution corresponds to a BTZ black hole with a short 
distance correction term, 




n(r) = \ -M + --1 , (15) 



and it describes a BTZ solution conformally coupled to a scalar field [22] . 

To introduce a brane we solve the corresponding junction conditions given 
by the boundary Einstein equations ((8| using the induced metric shown in 



(10). For the case when n{r) corresponds to the BTZ black hole (14), and 
the brane cosmological constant is given by A3 = — I /I 2 , we found that the 
energy-momentum tensor in (|8| is null. Therefore, the BTZ black hole is 



localized on the brane in vacuum. 



When n(r) is of the form given by (15), the energy momentum tensor 



necessary to sustain such a solution on the brane is given by 

2S = diag , (16) 



which is conserved on the brane [23] . 

These solutions extend the brane BTZ black hole into the bulk. The 
warp function f 2 (p) gives the shape of a 'throat' to the horizon, whose size 
is defined by the scale y/a, which is fine-tuned to the length scale of the 
five-dimensional AdS space. 



3 The Kerr-Schild Method 

In this section we will apply the generalized Kerr-Schild method to the met- 
rics analyzed in [13] to generate new solutions. This procedure consists on 
defining a new metric qmn starting from a known one c/mn as follows 

Qmn = Qmn + 2 H(r, p) t M i N , (17) 
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and checking if this new metric g fiu satisfies the Einstein equations. Here H 
is an arbitrary function of the coordinates r and p in this specific case, and 
£m is a null geodesic vector (in the background metric) described below (18). 



A. H. Taub [15] introduced this approach in the form just described. 
There is an extended family of solutions generated by this method. Among 
others, it reproduces the standard Einstein's solutions: Schwarzschild, Reissner- 
Nordstrom, and (obviously) the Kerr solutions starting from a flat, empty 
spacetime. We extend this formalism to the BTZ string-like solutions de- 
scribed in the previous section. 

Our case is quite different from those just described since it includes 
the Gauss-Bonnet term in the original action Q. The resulting equations 
are at most quadratic in the function H(r,p) due to the index contractions 
displayed in the GB term. 



3.1 Null Geodesic Vector 

Let us consider a vector £m = C-uij, p) which obeys the following conditions, 

0, £ M £ M = (18) 



where capital indices run over all the coordinates (t, r, <f>, p, 9). 

The definition of the following operator will show to be useful in solving 
these equations, 

d „, „ d 



D 



n 2 £ r ~ + /% 7t 
Or op 



(19) 



Using this operator and the geodesic equation ( |18| ), we obtain the follow- 

(20) 



ing equations for the t, <ft, and 9 components of the null vector £a 

1 



P 



D 



, for \ = t.d> and 9. 



The r and p components are more involved and read, 

n „o . „ Q 1 , 



7- Dt < + 7> 



— I ? + n h £l -i 

ti r 3 



j 2 H 



J 2 , 



n 2 C t U l r r 2 l 4> 



+ 



b 2 







0. 



(21) 
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In our convention over-dot means differentiating with respect to r while a 
prime corresponds to a derivative with respect to p. 

An additional constraint comes from £m being a null vector (18). This 
condition is 



1 



1)2 



rr 



+ 



2 , H 

' p b 2 



0. 



(22) 



The Kerr-Schild formalism requires to know at least one explicit solution 
for the null geodesic £a before we begin our search for a new solution of 
Einstein's equations. We start considering the following assumptions about 
the coordinate dependence of the radial brane and bulk components of the 
lightlike vector, £ r and £ p 



£ r (r), 



and 



(23) 



With these assumptions a solution for the set of Eqs.(|20j) is readily found. 
For the components i = t, 8, and <fi we have 



£i(r,p) = d exp Ktip 



dr 



n 2 £ r (r) 



x exp —Ki 



dp 



Pip) if 



(24) 



A simple non-trivial solution can be obtained setting K\ = 0. In this 
case the solutions are 



E, 



L , and 



K, 



(25) 



where E, L, and K are constants related to the energy and the brane and bulk 
components of the angular momentum of the particle following the geodesic. 

There are still two components left to solve £ r (r) and £ p (p). Introducing 
these expressions in the Eqs.(21 ) and doing the corresponding simplifications 



we arrive to 



n{rf {£D + 



E 2 + (n 2 )£l + 



L 1 



W)' + 



J 2 



\ E 2 + n 2 £l + \L 2 



ir 



0. (26) 
% = (27) 
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Let us work with the last equation. Employing the null vector condition 
(22) and factorizing conveniently we obtain an expression that can easily be 

K 2 ' 



integrated 



1 



f{pY 



f(pY W + 



b( P y 



o. 



from here we find the general solution for £ p (p), 

K 2 < 2 



Hp) 2 f( P Y 



(28) 



A new constant £ has been introduced here. If K ^ 0, Eq.(28) is valid 
as long as its right hand side remains positive and p > 0. 



We now solve t r from Eq.(26). This equation can be written as a total 



derivative and after one integration becomes 



E 1 



+ n(r) 2 £ 2 r + 



L 2 



X' 



where \ is a constant of integration. However, this constant is not a new 



parameter since this component must fit in the null vector restriction Eq.( 22 ) 
Inserting this expression in Eq.(22) we obtain that x 2 + £ 2 = 0, fixing \- 



With this last step we have solved Eqs.(18) for each of the components 
of the null geodesic vector. 



The null geodesic vector takes the following general expression, 



'M 



E 



n- 



E 2 — 



' L 2 



p 



K 2 



K 



(29) 



3.2 Solutions Generated by the Kerr-Schild Method 

We choose the following null geodesic vector, 



n- 



^n 2 ,L,0,0 



(30) 



where L is a constant associated to the angular momentum of the test parti- 
cle, and a function H of the form H(r, p) = h\{r) h^p)- We use this Ansatz 



in (17) and replace back into the Einstein equations obtaining eight non-null 
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equations. Factorizing (pt) and (pr) components we arrive to the following 
equation 

[2h 2 (p)f-fh' 2 (p))[h 1 (r)+h , l (r)r) = 0. (31) 



Once we solve Eq.(31 ), the rest of the equations are automatically fulfilled 



provided that the corresponding constraints in Table [T] are satisfied. 

If we choose to solve this equation for hi(r), we obtain hi(r) — a/r, where 



a is a constant, and h 2 (r) remains free. When putting back into (17) the 
new metric turns into 



ds 2 



2„2 , o u \ dt 2 _ 4 2_ dUr + IJ +2 1 )dr 2 

\n 2 r n 4 I 



-f'n z + 2-h 2 
+f 2 r 2 d<j)' z + dp' + tfd6 



r / r n 

2 , J„2 , >2jfl2 



(32) 



This metric represents a 5-dimensional gravity solution, which can be diag- 
onalized in few cases depending on the form of h 2 (p). Notice that the line 



element (32) does not include a boundary membrane. As our main interest 



here is to find braneworld solutions, naturally our next step is to embed a 
brane in this bulk metric. In order to proceed to this point, the first require- 
ment we find is that the induced metric must fulfill the junction conditions, 
i.e., the 3-dimensional Einstein equations ([8]). As far as we know, the BTZ 
metric is the only solution describing a black hole in a (2+1) spacetime, thus, 
h 2 (p) becomes constrained to a multiple of the warp factor f 2 (p) to be able 
to recover a BTZ-like metric on the brane. In particular, if h 2 (p) = f 2 {p), 
we can make a coordinate transformation to end up with the BTZ string 
coupled to a brane scalar field as it will be shown below. 

Alternatively, if we solve Eq.(31) for h 2 (p), we get h 2 (p) = f 2 (p), and 
hi(r) becomes arbitrary. In this case the choices for hi{r) are several and 
give place to the following new solutions. Nevertheless, as the presence of 
the brane imposes junction conditions we should point out that the 
right options for hi(r) will be determined by the requirement of yielding a 
physically meaningful brane energy-momentum tensor. 



3.2.1 Charged BTZ string 

Let us first consider the solution f(p) = cosh(p/2-^/a) and b(p) = 2/3-y/a smh.(p/2^/a). 
By choosing h\(r) = ^- In r and L = we arrive to the following metric, 



ds 2 = f 



-(n 2 -Q 2 \n r)dt 2 + 



2Q 2 In r 



rr 



dt dr + 



n 2 + Q 2 ln r s 



dr + r 
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+dp 2 + b 2 d6 2 . (33) 

In order to obtain a more familiar form of the metric we make the following 
coordinate transformation, 

Q 2 In r dr . . 

dt = dt + 2 ( 2 n21 \ • 34 

n z {n z — Q 2 m r) 



This change cancels out the non-diagonal term such that the metric (33) 
turns into 

ds 2 = f 2 (-n 2 dt* + ^ + r'd^ + dp 2 + b 2 d6 2 , (35) 

with n 2 = —M + r 2 /l 2 — Q 2 In r, which describes a charged BTZ string whose 
charge is confined to the brane. In order to verify this statement we calculate 
the brane energy momentum tensor necessary to hold this solution and we 
obtain 

This is precisely the stress energy tensor related to a charged object in (2+1) 
dimensions. 



3.2.2 BTZ string coupled to a brane scalar field 

Working with the same expressions for f(p) and b(p) we now choose h\{r) 
S- and the resulting metric is 



ds 2 



— \n — 



dt 2 + -^-dt dr+\ 



n 2 r 



n A r 2 



+dp 2 + b 2 d9 2 . 
We now make the coordinate transformation 

c 



dt = dt + 



n 2 r(n 2 — C/ r ) 



to arrive to the following metric 
ds 2 = f 2 (-n 2 dt 2 + 



dr 2 



n- 



dr + r 



2.i2 



+ dp 2 + b 2 d6 2 



(37) 



(38) 



(39) 



where n 2 = —M+r 2 /l 2 — (/r, which corresponds to a BTZ black hole coupled 
to a scalar field on the brane. This solution had already been found in [13]. 
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3.2.3 Charged BTZ string coupled to a brane scalar field 



Another possible combination is to choose the original metric to be (35) and 
hi{r) 

-M + r 2 /l 2 



C/2r. Following the same procedure as the previous case, we find 
the same metric as (39), but with n 2 



Q In r — Qjr. If we 



compute the energy-momentum tensor on the brane, we find 

Q 2 O 



diag 



o 2 _c_ 

2r 2 2r 3 ' 



2r 2 2r 3 



2r 2 



(40) 



which corresponds to a charged BTZ black hole coupled to a scalar field on 
the brane. 



3.2.4 BTZ string with angular momentum 



In order to add angular momentum to the original BTZ string, we pick 
hi{r) = c, where c is a constant. In this case L ^ 0. Thus, the metric takes 
the following form, 



— in 



2c) dt 2 + 



4c 



L 5 



n 2 dt dr + AcL dt < 



n- 



+ 



n 



2 r 2 + 2cr 2 - 2cn 2 L 2 



n^r 2 



-dr 2 + 



AcL 



L< 



n L 



n 2 dr dcf) 



+ (r 2 + 2cL 2 )d(p 2 ] + dp 2 + b 2 d6 2 . 

Introducing the following transformations, 

dt = di + u(r)dr 
d(p = d(p + v(r)dr , 



with 



u(r) 
v{r) 



2cr 2 



L 2 



n 2 (n 2 r 2 — 2r 2 c + 2cL 2 n 2 ) 



n 



L—u(r) . 



(41) 

(42) 
(43) 

(44) 
(45) 



Moreover, if we define R 2 = r 2 + 2cL 2 , J = -AL, and M = M + 2c(L 2 /I 2 + 1), 
we arrive to the metric for a rotating BTZ-string, 



ds 2 = f 



-tfdP + ^ + R 2 



-J 



i^2B 2 



dt 



+ dp z + b 2 d6 2 



(46) 
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where n 2 (r) = —M + R 2 /l 2 + J 2 /AR 2 . The corresponding energy- momentum 
tensor on the brane can be calculated, and we find that it vanishes, as it 
should be for a rotating BTZ brane black hole. 

Analogously, when we use the solutions f(p) = 1, b(p) = 7sinh(p/7), 
and f(p) = 1, b(p) = 2(3jasmh(p/2Ja), we also arrive to several solutions 
involving angular momentum and scalar fields. We should stress that in these 
cases we did not find any charged solution. 



4 Geodesic Structure 

In this section we study the geodesic behaviour in the background of the 
solutions displayed in Table [1} 

Let us begin our study by considering the Lagrangian for the BTZ black 
hole in codimension-2 branes with the solution f(p) = cosh (§3^) and b(p) = 

2/VSsirm(4 s ), 



C = cosh 



n(r)V + 



+ r 2 4> 2 I + p z + 4/3' z a sinh 



n(r) 2 



2^. 



(47) 

where a dot indicates derivative with respect to the affine parameter A. 

As it is independent of t, 6 and 0, we can write the following equations 
of motion, 

\ 2 

p > n(r) 2 i = -2E, 



DC 

dC 
d$ 
dC 

~de 



—2 cosh 
2 cosh 



\ 2 
P \ 



2 Jet 



8(3 2 a sinh 



P 



2J^ / 



r 2 <f) = 2L, 
\ 2 

e = 2K. 



(48) 
(49) 
(50) 



where L and K are the angular momenta of the particle related to and 
9 coordinates, respectively. Notice that Eq.(48) gives us a relation between 
the coordinate time t and the affine parameter A. 

With these equations the Lagrangian (47) can be written as 

—E 2 r 



cosh 



2J^ / 



+ 



cosh 4 (^) n{rf n(r) 2 r 2 cosh 4 (^) 
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where h = or —1 is a parameter describing lightlike and timelike geodesies, 
respectively. 



4.1 Geodesies on the Brane 



Let us consider a particle with K = 0. Using Eq.(51) at the position of the 
brane (p = 0) we can find the effective potential for the geodesic motion, 

r 2 = £ W(r) - />) =► V: ff = n\r)(^- h y (52) 

This equation can be integrated to obtain the orbits as well, 



dr 



E 2 — n 2 {r) (~ — h\ . (53) 



4.1.1 BTZ Case 

For radial geodesies (L = 0) the effective potential becomes 



y»2 



V e ) f = n\r)h=\-M + -jh. (54) 

In this case lightlike geodesies are just straight lines. 

For the timelike case, we display the potential in FigjTJ The corresponding 
orbits can be obtained by direct integration of Eq.([53|), 



r(A)=/^^[l + sin(2A)] 1/2 , (55) 

and they are shown in FigjTJ We see that geodesies have an oscillatory 
behaviour outside the event horizon, however, there are no stable orbits since 
some of the particles could cross it and never return. Notice that particles 
with E < are not allowed since V e fj = at the event horizon. 

For particles with angular momentum the effective potential turns out to 

be 

V eff = i-M + T pj(^-h). (56) 
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Figure 1: Effective potential (left) and orbits (right) for radial timelike par- 
ticles on the brane. These graphs correspond to I — 1, M — 0, 5. Notice that 
geodesies have an oscillatory behaviour outside the event horizon, however, 
there are no stable orbits since some of the particles could cross it and never 
return. Particles with E < are not allowed since V e ff = at the event 
horizon. 



We plot this potential for lightlike and timelike cases in Fig|2j where we 
set M — 0, 5, I — 1, L — 2, and L = 5. In both cases the potentials intercept 
themselves at the event horizon. Remark that far from the event horizon 
the lightlike potential has an asymptotic behaviour that can be inferred from 
Eq.Q, i.e., V eff ->• L/l. 

In order to obtain the orbits for the lightlike case (h = 0) with angular 



momentum we can integrate Eq. ( 53 ) to attain 



r(A) = ± 



\ 



EP-^A A 2 ML2 

P J E 2 — L 2 /l 2 



(57) 



and we should stress that only the plus sign has physical meaning. 



Analogously, we perform the integration of Eq.(53) in the timelike case 
and we arrive to 

1/2 



r(A) 



IE 2 l 2 + Ml 2 - L 2 



1 



ANIL 2 



\ (E 2 + M-L 2 /l 2 ) 2 l- 



sin(2A) 



(58) 

The corresponding orbits are presented in Fig|3] From this figure we can 
see that lightlike geodesies with low energies fall unavoidably into the event 
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Figure 2: Effective potential for lightlike (left) and timelike (right) brane 
particles with different angular momenta (L = 1, 2, 5). Remark that all the 
potentials cross themselves at the event horizon. The lightlike potential has 
an asymptote given by L/l, while the timelike potential grows with no limit. 



horizon. However, particles with high energies can escape from the black hole 
showing that energy extraction is possible but only with massless particles. 
On the other side, the timelike orbits have basically the same shape as those 
in the radial case (Fig(l|D), the only effect of L is to increase the amplitude 
of the oscillation. Again some of them can cross the event horizon depending 
on the energy of the oscillation. In both lightlike and timelike similar 
qualitative behaviour was also found in the study of pure BTZ geodesic 
structure [21]. 



4.1.2 BTZ with Electric Charge and Scalar Field 

In this case the resulting potential when adding charge or a scalar field to 
the BTZ solution is again given by Eq.(52) with the corresponding n(r) as 
follows, 

n(r) = W -M + ^ - Q 2 ln(r) , charged BTZ (59) 



n(r) = \I-M + — --, BTZ + scalar field (60) 



/ 

Both potentials are shown in Fig|4j 
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Figure 3: Orbits for lightlike (left) and timelike (right) brane geodesies with 
angular momentum (L = 2). The energies of the particles are shown in the 
legend. Notice that in the lightlike case orbits with low energy fall into the 
event horizon whereas particles with high energy can escape from the black 
hole. For timelike geodesies the orbits have basically the same shape as those 
in the radial case (Figflb), the only effect of L is to increase the amplitude 
of the oscillation. 




Figure 4: Effective potential for timelike geodesies on the brane for BTZ 
with charge (left) and with scalar field (right). These graphs correspond to 
Q = 1, 2, 3 and £ = 1,5, 10. In the former case the potential admits oscillating 
orbits that can fall into the event horizon when Q is small, however, as Q 
grows, the minimum of the potential is shifted outside the horizon making 
possible the existence of stable oscillating or bounded geodesies. In the latter 
case only unstable oscillations are allowed. 
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We can deduce from these graphics that the charged BTZ potential ad- 
mits oscillating orbits that can fall into the event horizon when Q is small, 
however, as Q grows, we see that the minimum of the potential is shifted 
outside the horizon making possible the existence of stable oscillating or 
bounded geodesies. The potential corresponding to BTZ coupled to a scalar 
field shows that just unstable oscillating orbits are allowed. 



4.2 Geodesies in the Bulk 

Here we study the geodesies that explore the extra dimensions. Although 
standard particles are not allowed to travel outside the brane, we perform 
this analysis as a way to acquire a better understanding of the geometry 
of these solutions. For this analysis we will consider the r coordinate lying 
outside the black hole horizon, for instance at r = 2yMl = 2r#. In this case 
it will be convenient to write Eq.(51) in a different way, 

E 2 



P 



h + 



L'- 



K 2 



Mcosh 2 (^) 



2MPcosh 2 (^) 



4/3% sinh 2 (^) 



Defining a new variable u as 



(61) 



u 



2 yfa. sinh 



a 



and replacing in Eq.(61) it becomes 

L 2 



u 



e 2 - 



2Mf 



+ 1 + 



w 2 \ / K 2 

4a \{3 2 u 2 



h 



where e 2 = E 2 /M. Thus, we can define an effective potential given by 



2M1 2 



+ 1 



u 
4a 



' K 2 
~Wu 2 



h 



(62) 



(63) 



(64) 



For radial geodesies (L = K = 0) notice that as the effective potential 
vanishes in the lightlike case, the orbits are just straight lines. On the other 
side, the orbits for timelike geodesies can be found by integrating Eq.(63) 
and replacing u from Eq.(62). Thus, we obtain, 

V E 2 — 1 sin 



p(A) = 2y/a arcsinh 



A 



a 



(65) 
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Some of the orbits are depicted in Fig|5j This graph shows oscillating tra- 
jectories that cross the brane. This implies that particles leaving the brane 
can return in the future. This fact opens up the possibility for the existence 
of shortcuts, paths connecting two points which are shorter in the bulk than 
on the brane [23] ■ 111 addition, note that the trajectory with the energy that 
corresponds to the minimum of the effective potential is on the brane. 




Figure 5: Effective potential and timelike orbits for radial particles in the 
bulk. From the oscillating behaviour of the trajectories we can infer the 
existence of shortcuts since particles leaving the brane can return to it. Note 
that the trajectory with the energy that corresponds to the minimum of the 
effective potential (E = 1) is entirely on the brane. 

Now we turn to the case of timelike and lightlike geodesies with K ^ 
and/or L ^ 0. 

In the lightlike case, if K = and L ^ we have a constant potential 
V e ff = 2M2- The corresponding timelike case has the same potential as the 
L = geodesies, only shifted by a constant, thus, the shape of the geodesies 
is the same as the one shown in Fig|5]D. 

When K ^ and L = 0, we obtain the orbits displayed in FigJH} We 
can notice that in the lightlike case the particles seem to be scattered by 
a barrier-like potential near the brane. Regarding the timelike geodesies, 
we can see an oscillating behaviour around certain position parallel to the 
brane, which corresponds to the orbit with the minimal permitted energy. 
Additionally, the higher the energies of the particles are, the closer to the 
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brane they can reach, but never cross it. 
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Figure 6: Orbits for lightlike (left) and timelike (right) bulk particles with 
angular momentum (K = 2, L = 0). Lightlike geodesies show a scattering 
behaviour. While timelike particles oscillate near the brane, and the higher 
their energies are, the closer to the brane they can reach, but never cross it. 



Forthwith, let us check the geodesic behaviour for the solution /(p) = 1 
and b(p) = 7sinh(p/7) (see Table [T]). By fixing r we can write an analogous 
equation to (|5lj), 



E 2 L 2 . 2 K 2 



C = — k~t\ + ~r + P + 9 . ,2/ / x = h ■ ( 66 ) 

n z (r) r z j A smh (p/7) 



So that, 



E 2 L 2 K 2 



P = "TTT V 9-i2/ / \ + h ■ (67) 

n z (r) r l 7 2 smh (p/7) 

Fixing r the effective potential becomes, 

L 2 K 2 
V "ff = ^MP + 72 sinh 2 (p/ 7 ) " k " (68) 

Notice that this potential is constant when K = 0. The corresponding graph 
can be seen in Fig|7j Observe that timelike and lightlike potentials mark off 
from each other just by a constant that shifts the entire curve. 
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Figure 7: Bulk effective potential for / = 1, with L = 0, 2 and K = 1. We 
can see that close to the brane the potential has a scattering behaviour, while 
far from the brane it becomes constant, so that any bulk geodesic behaves 
like a free particle. 

This is a scattering potential where no other behaviour is possible because 
of the infinite asymptotic barrier caused by the diverging term. We can infer 
that gravitational signals or particles already in the bulk cannot reach the 
brane and those originated on the brane cannot travel far from it. 

4.3 BTZ with Angular Momentum 

In this subsection we study the particular case of geodesies on the brane and 
in the bulk for the BTZ black hole with angular momentum in co-dimension 
2-brane. The Lagrangian can be written as follows, 



C = f 



-AM + -nr + -z\1? + 



4 I p r 2 J (-4M + 4r 2 /P + J 2 /r 2 ) 




+ p 2 + b 2 e\ (69) 



and the equations of motion become 
dC „ „ 2 2 , Ui 



2f 2 r 2 (--- + ( j ) ) = 2L, (70) 
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dC 
dC 



p 



1 ( 4r2 

-- -4M + — + — 



; 2 



t - 



1 Jt 

2r2 



J 



26 2 fl = 2K . 



-2E, (71) 
(72) 



With these equations the Lagrangian becomes 
1 E 2 — LJE/r 2 + ML 2 /r 2 — L 2 /l 2 



h 



P 



-M + r 2 /l 2 + J 2 /Ar 2 



Pr 2 . 2 K 2 

-M + r 2 /l 2 + J 2 /4r 2 6 2 ' 

(73) 

First, we look for the geodesies on the brane. Solving this equation for r we 
obtain, 



^ r 2 J 2 

J 2 \ P + 4T 2 



K 2 ^ 
b 2 



+ p 2 -h 



+ 



.Z 1 



E — 



LJE ML 2 



L 2 ' 
/ 2 



(74) 



In this case we will go through the integration to find the orbits directly. As 
the brane is located at p = 0, /(0) = 1 in any solution. In addition, we set 
K = to avoid singularities. Thus, the integral for the orbit becomes, 

dr 



A-A n 

-M + r 2 /l 2 + J 2 /4r 2 )h + (E 2 - LJE/r 2 + ML 2 /r 2 - L 2 /l 2 ) 

(75) 

For the timelike case (h — — 1) and setting Aq = the orbits are given by 



r(A) 



± 



sin(2A//) 



I 2 ( 2 L 2 

+ 2 [ M + E ~P 



{Ml 2 + L 2 + E 2 l 2 ) 2 - {2 ELI + U) 2 

1/2 



(76) 



For the lightlike case (h = 0) we have 



r(A) = ± 



E 1 



T 2 



(ML 2 - LJE) 
E 2 -L 2 /l 2 



1/2 



(77) 



The possible orbits are displayed in FigjSJ These graphs show that geodesies 
can cross the event horizon in both lightlike and timelike cases. According 
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to the energy of the particle, it can fall directly or after a roundabout in the 
region exterior to the horizon. Remark that the circular orbits corresponding 
to E = wind inside the ergosphere while the other orbits with E > extend 
farther. In addition, massless particles with high energies can escape from 
the gravitational attraction of the black hole allowing energy extraction, a 
behaviour already noticed in the pure (2 + 1) BTZ case [24]. 




Figure 8: Lightlike (left) and timelike (right) brane geodesies for BTZ so- 
lution with angular momentum. Both positive and negative branches are 
plotted. The energies of each geodesic are displayed on the legend. No- 
tice that in both cases geodesies can cross the horizon, some of them after 
a roundabout and some others directly. Only in the lightlike case particles 
with high energies can escape from the black hole allowing energy extraction. 



Now we turn to the geodesic motion in the bulk. Accordingly, we write 

Eq.(73) taking r = and r = \J2Ml 2 + 2ly/M 2 l 2 - J 2 > r E , where r E is 
the ergosphere position. We first check the case f(p) = cosh(p/2-y/a), and 
b(p) = Ifi^fa sinh(p/2-y/a). After integrating, the orbits turn out to be 



P(A) J 



2^/a arcsinh 



2^/a arcsinh 




4C 2 a + \ 2 C{ 



C\CL 



Ci + VCT - 4C 2 sin 



A 



a 



(78) 



(79) 



for lightlike and timelike geodesies, respectively. The constants C\ and C2 
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are given by 



C 1 = h-C 2 + \ 



llsjMH 2 - J 2 (E 2 - L 2 /l 2 ) + 2l 2 E 2 M - L 2 M - JLE~ 

8Ml(Ml + y/M 2 l 2 - J 2 ) - 5 J 2 

(80) 



Some trajectories are shown in Figj9] (upper graphs). In both cases although 
the geodesies approach the brane, they cannot cross it. Analogously to the 
non-rotating case, we see that the brane acts like a repulsive barrier that 
blocks the exchange of signals between brane and bulk. Nevertheless, whereas 
lightlike geodesies are scattered when trying to reach the brane, timelike 
trajectories oscillate about a fix distance from the brane. 

Similarly, the trajectories corresponding to the solutions f(p) = 1, b(p) = 
7sinh(p/7) and f(p) = 1, b(p) = 2(3y/a smh(p/2y/a) undergo an analogous 
pattern displayed in Fig|9] (lower graphs). Again we see a repulsive barrier 
behaviour at the position of the brane that scatters lightlike and timelike 
particles. This prevents both signals originated on the brane to permeate 
the bulk and signals already in the bulk to arrive to the brane. 



5 Conclusions 

In this paper we considered the black hole solutions of five dimensional grav- 
ity with a Gauss-Bonnet term in the bulk and an induced gravity term on a 
2-brane of codimension-2 [T3]. In order to explore the possibility of having 
more solutions to the Einstein equations we applied the Kerr-Shild method 
using as a background the above-mentioned solutions. This method adds to 
the original metric a term depending on a null geodesic vector and a scalar 
function H . Due to the presence of the GB term the resulting equations 
are at most quadratic in H, but still solvable by choosing the appropriate 
geodesic vector. Working out the Einstein equations for the function H we 
found additional solutions, which include charge, angular momentum, and 
brane scalar fields coupled to the brane black hole. The corresponding brane 
energy-momentum tensors were also computed. These results lead us to 
deduce that the solution with free n(r) found in [13] does not have any 
constraint of diagonality. 
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Figure 9: Lightlike (left) and timelike (right) geodesies in the bulk for BTZ 
with angular momentum case with / = cosh(p/2\fa) (up) and / = 1 (down). 
In both cases we observe that particles cannot reach the brane. We can infer 
that the brane acts as a repulsive barrier both confining particles already on 
the brane and preventing the entrance of bulk particles. 

Furthermore, we studied the geodesic behaviour in the background of the 
original and new solutions. For this purpose we solved the Euler-Lagrange 
equations for the variational problem associated with the metric. Among 
our results we can distinguish different cases according to the energy of the 
particle, the associated effective potential, and the components of the angular 
momentum. 

In the case of geodesies on the brane we can discriminate two main cases, 
timelike and lightlike geodesies. The timelike orbits, in general, display os- 
cillations that can become unstable making a particle cross the event hori- 
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zon and never return. Nevertheless, the charged BTZ string-like solution 
shows an additional behaviour. When the charge is small, the paths are the 
same as those appearing in uncharged or scalar field coupled solutions; how- 
ever, as the charge grows, the minimum of the effective potential is shifted 
outside the event horizon and makes possible the existence of stable oscil- 
lations or bounded orbits for particles with low or negative energy. The 
lightlike geodesies exhibit particles with low energies crossing the event hori- 
zon. Nonetheless, particles with high energies can escape from the black hole 
and, thus, allow the extraction of energy, a result already noticed in pure 
(2+1) BTZ black holes geodesic structure [24J. 

Concerning the geodesies in the bulk, we fixed r > and studied two 
main cases, namely, / = cosh(p/2y / a) and f — 1. In the first case, we ob- 
served that timelike trajectories corresponding to particles without angular 
momentum are oscillating about the position of the brane. This implies that 
a particle leaving the brane can return to it opening up the possibility of 
shortcuts, i.e., paths that are shorter in the bulk than on the brane [25J. 
The same behaviour is verified for particles having an angular momentum 
component along the brane (L 7^ 0). When K (the particle's angular momen- 
tum component along the bulk) is switched on, lightlike geodesies encounter 
a barrier-like potential around the brane. Moreover, although timelike tra- 
jectories are still oscillating and approach the brane more and more as their 
energy increases, they never cross it. This result prevents the exchange of 
bulk and brane signals, i.e., particles living in the bulk do not enter the 
brane, and particles residing on the brane cannot go far into the bulk. In the 
case / = 1, both timelike and lightlike particles undergo scattering near the 
brane, but move in straight lines far from it because the effective potential 
becomes constant at this far region. 

The case of a BTZ black string-like object with angular momentum was 
treated separately due to the non-diagonal terms appearing in the metric. We 
found that both lightlike and timelike geodesies on the brane can cross the 
horizon directly or after a roundabout in the ergosphere (particles with the 
lowest allowed energy) or even outer regions. However, massless particles 
with high energies are able to escape allowing energy extraction from the 
black hole. With respect to bulk geodesies, we determined that regardless 
the expression for /, the effective potential has a barrier-like structure around 
the brane, which both yields scattering bulk geodesies and confines particles 
already on the brane. In particular, timelike geodesies in the background of 
/ = cos\i{p/2yfa) oscillate around a fixed distance from the brane but never 
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make contact with it. 

It would be interesting to investigate if these features remain valid in 
codimension-2 (3 + 1) brane scenarios, and if these results may help uncover 
other interactions between the brane and the bulk. However, as this is out 
of the scope of this paper, we expect to address these questions in a future 
work. 

Acknowledgments 

B. C-M. and S. A. thank the hospitality of Facultad de Ciencias Fisicas y 
Matematicas of Universidad de Chile, where this work was accomplished. 

References 

[1] A. Chamblin, S. W. Hawking and H. S. Reall, Phys. Rev. D 61, 065007 
(2000) 

[2] R. Gregory, Class. Quant. Grav. 17, L125 (2000); R. Gregory and 
R. Laflamme, Phys. Rev. Lett. 70, 2837 (1993); G. Gibbons and 
S. A. Hartnoll, Phys. Rev. D 66, 064024 (2002). 

[3] T. Shiromizu and M. Shibata, Phys. Rev. D 62, 127502 (2000); 
A. Chamblin, H. S. Reall, H. A. Shinkai and T. Shiromizu, Phys. Rev. 
D 63, 064015 (2001); T. Wiseman, Phys. Rev. D 65, 124007 (2002). 

[4] T. Shiromizu, K. I. Maeda and M. Sasaki, Phys. Rev. D 62, 024012 
(2000); A. N. Aliev and A. E. Gumrukcuoglu, Class. Quant. Grav. 21, 
5081 (2004). 

[5] R. Casadio, A. Fabbri and L. Mazzacurati, Phys. Rev. D 65, 084040 
(2002); K. A. Bronnikov, V. N. Melnikov and H. Dehnen, Phys. Rev. D 
68, 024025 (2003). P. Kanti and K. Tamvakis, Phys. Rev. D 65, 084010 
(2002). 

[6] N. Dadhich, R. Maartens, P. Papadopoulos and V. Rezania, Phys. Lett. 
B 487, 1 (2000); M. Bruni, C. Germani and R. Maartens, Phys. Rev. 
Lett. 87, 231302 (2001). G. Kofinas, E. Papantonopoulos and I. Pappa, 
Phys. Rev. D 66, 104014 (2002); G. Kofinas, E. Papantonopoulos 



28 



and V. Zamarias, Phys. Rev. D 66, 104028 (2002); A. N. Aliev and 
A. E. Gumrukcuoglu, Phys. Rev. D 71, 104027 (2005). 

[7] M. Banados, C. Teitelboim, and J. Zanelli, Phys. Rev. Lett. 69, 1849 
(1992). 

[8] R. Emparan, G. T. Horowitz and R. C. Myers, JHEP 0001, 007 (2000); 
R. Emparan, G. T. Horowitz and R. C. Myers, JHEP 0001, 021 (2000). 

[9] N. Kaloper and D. Kiley, JHEP 0603, 077 (2006). 

[10] M. Aryal, L. H. Ford and A. Vilenkin, Phys. Rev. D 34, 2263 (1986). 

[11] D. Kiley, Phys. Rev. D 76, 126002 (2007). 

[12] D. Dai, N. Kaloper, G.D. Starkman, D. Stojkovic Phys. Rev. D 75, 
024043 (2007). 

[13] B. Cuadros-Melgar, E. Papantonopoulos, M. Tsoukalas, and V. Za- 
marias, Phys. Rev. Lett. 100, 221601 (2008). 

[14] B. Cuadros-Melgar, E. Papantonopoulos, M. Tsoukalas, and V. Za- 
marias, Nucl. Phys. B 810, 246 (2009). 

[15] A. H. Taub, Ann. of Phys. 220, 326 (1981). 

[16] S. S. Komissarov, Mon. Not. R. Astron. Soc. 326, L41 (2001). 

[17] M. F. Huq, M. W. Choptuik, and R. Matzner, Phys. Rev. D 66, 084024 
(2002). 

[18] R. Alonso, N. Zamorano, Phys. Rev. D 35, 1798 (1987); L. Pena, 
N. Zamorano, internal publication (2005). 

[19] A. Anabalon, N. Deruelle, Y. Morisawa, J. Oliva, M. Sasaki, D. Tempo, 
R. Troncoso, Class. Quant. Grav. 26, 065002 (2009). 

[20] J. Ponce de Leon, Int. J. Mod. Phys. D 12, 757 (2003); W. Mueck, 
K. S. Viswanathan, and I. V. Volovich, Phys. Rev. D 62, 105019 (2000); 
D. Youm, Mod. Phys. Lett. A 16, 2371 (2001); F. Dahia, C. Romero, 
L. F. P. Silva, and R. Tavakol, J. Math. Phys. 48, 072501 (2007); S. Das, 
S. Ghosh, J-W. van Holten, and S. Pal, JHEP 0904, 115 (2009). 



29 



[21] P. Bostock, R. Gregory, I. Navarro, and J. Santiago, Phys. Rev. Lett. 
92, 221601 (2004). 

[22] C. Martinez and J. Zanelli, Phys. Rev. D 54, 3830 (1996). 

[23] G. Kofinas, Class. Quant. Grav. 22, L47 (2005). 

[24] N. Cruz, C. Martinez, and L. Peiia, Class. Quant. Grav. 11, 2731 (1994). 

[25] H. Ishihara, Phys. Rev. Lett. 86, 381 (2001); R. Caldwell and D. Lan- 
glois, Phys. Lett. B 511, 129 (2001); E. Abdalla, A. G. Casali, 
B. Cuadros-Melgar, Int. J. Theor. Phys. 43, 801 (2004). 



30 



